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We describe a method for measuring CP-violation parameters from which the Cabibbo-Kobayashi- 
Maskawa angle 7 may be extracted. The method makes use of the total decay rates in DK^ 
decays, where the neutral D meson decays to multibody final states. We analyze the error of the 
method using experimental CP-violation analysis variables that enable straightforward sensitivity 
comparison with other methods for extracting 7, and discuss the use of _B-factory and charm-factory 
data to obtain the relevant charm decay information needed for this measurement. Measurement 
sensitivities are estimated for the currently available B-factory data sample, and D decay modes for 
which use of this method can make a significant contribution toward reducing the total error on 7 
are identified. 
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I. INTRODUCTION 

An important part of the program to study CP vi- 
olation is the measurement of the angle 7 = 
^^E{~^ud^ub/^cd^cb) of the unitarity triangle related to 
the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing 
matrix [1]. Measurement of 7 performed with tree-level 
processes defines an experimentally allowed region for the 
apex of the unitarity triangle. This region should overlap 
with the region obtained from B° - B° and Bs - Bs mix- 
ing, assuming there are no significant new-physics contri- 
butions in the mixing amplitude. With this assumption, 
current Tevatron measurements [2] of the B^ — Bs mixing 
rate yield an indirect constraint on 7 that is much tighter 
than direct measurements [3]. Therefore, precise direct 
determination of 7 presents an opportunity to conduct 
an accurate test of the Standard Model. 

The decays B — > DK can be used to measure 7 with 
essentially no hadronic uncertainties, exploiting interfer- 
ence between the b ucs and b —^ cus amplitudes of 
the decays B D^K and B D^K, respectively [4]. 
Interference takes place when the D meson^ is observed 
in a final state F that is accessible to both and 
decays. Such measurements can be conducted with quite 
a few D and B decay modes, including those with excited 
charm and strange mesons, involving different methods 
for constructing and optimizing CP- violation observables 
and measuring parameters related to 7. In fact, there has 
been a healthy stream of new ideas in this area since the 
basic method was first proposed in 1991 [4]. The dif- 
ferent parameters of the various measurements are then 
combined statistically, yielding confidence intervals for 
7 [3]. The statistical sensitivity provided by each mode 
and method is generally poor, mainly due to the strong 



^ We use the symbol D to indicate any linear combination of a D' 
and a meson state. 



CKM suppression (and, for most modes [5], color sup- 
pression) in the b ucs transition. This necessitates the 
exploitation of as many modes and methods as possible, 
in order to achieve a small combined error on 7. 

The most accurate 7 measurement method to date de- 
termines 7 by analyzing the Z?-decay event distribution 
in B± DK^ with multibody D decays [6, 7]. This 
method was initially applied to the Cabibbo-favored de- 
cay D K°TT+Tr^ [8, 9], and the BaBAR Collabora- 
tion later used it with K'^K'^K~ [10] and the Cabibbo- 
suppressed decay D tt+tt'tt" [11]. A simulation 
study has also been conducted for the four-body mode 
D K+K-TT+n- [12]. 

As originally proposed [6] , this method extracts the an- 
gle 7 from measurements of dT^{P)/dP, the differential 
decay rates of B^ DK^ at each phase-space point P 
of the multibody Z?-decay final state F. However, mea- 
surements done with the final states F = if^Tr+Tr" and 
F = K'^K~^K~ have only made use of the phase-space 
distributions, given by the relative differential rates 

dt^jp) _ dr^jp) 1 
dP ~ dP r^' ^' 

where 

p f dTl(P) 
^±^J -^T^dP (2) 

are the total decay rates. Thus, these measurements 
were sensitive only to the dependence of the rates on the 
point P, not to their integrated values F^. By contrast, 
the BABAR measurement with F = 7r^7r~7r° used both 
dr^{P)/dP and F^. For that mode, the total decay 
rate F^ gave more precise information about the CP- 
violation parameters than the phase-space distribution 
dT^{P)/dP. While most measurements and sensitivity 
estimates have focused on use of dT^{P) / dP for learn- 
ing about 7, it is important to identify and study the 
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decay modes for which the total decay rate has competi- 
tive sensitivity to the CP-violation parameters. This will 
help ensure that all useful modes are utilized for measur- 
ing 7, while preventing much effort from being wasted on 
data analysis of decay modes that are not promising. 

The purpose of this paper is to provide the tools for es- 
timating the CP-parameter sensitivities of measurements 
of the absolute decay rates for different D decay 
modes. We demonstrate that a good estimate of the sen- 
sitivities is provided by a single mode-dependent param- 
eter. The impact of each mode on the combined error of 
7 depends on values of strong phases and decay distribu- 
tions that in many cases are not well known yet. How- 
ever, our general analysis of the sensitivities, performed 
in terms of CP-violation parameters similar to those used 
in the most accurate experimental analyses to date, pro- 
vides a good indication as to when using the integrated 
decay rates is expected to improve the overall precision 
on 7. Since the combined error on 7 depends on many 
measurements, its full estimation is not within the scope 
of this paper and is not attempted here. Rather, we com- 
pare the sensitivity of the absolute-rate analysis to that 
of the current-best phasc-space-distribution analysis us- 
ing comparable experimental CP-violation variables. 

We present the formalism for the decay rates mB~ ^ 
DK~ with multibody D decays in Section II. Meth- 
ods for measuring important charm-decay quantities are 
discussed in Section III. The sensitivities with which 
the CP-violation parameters are obtained from the to- 
tal rates are calculated in Section IV, then estimated for 
self-conjugate D decay modes in Section IV A and for 
non-sclf-conjugate modes in Section IV B. We provide 
numerical estimates for several cases, in which enough in- 
formation is available for carrying out this calculation, in- 
dicating the promising and not-so-promising final states 
for this type of analysis. Actual data analysis of the type 
discussed here has been performed for only one of the de- 
cay modes we study, D 7r+7r~7r°. For all other modes, 
the estimates we provide are new. 



II. 



DK- DECAY RATES 



Consider the decay DK^_^ D F{P), where D 

is a superposition of the and states, F represents 
the particles comprising a multibody final state accessible 
through both D'^ and decays, and P is a specific point 
in the phase space of F . We are also interested in events 
involving the decay D F{P), where F{P) is the CP 
conjugate of F{P). The B-meson decay amplitudes to 
final states with specific charm flavor are parameterized 
and denoted in this paper in the following way: 

A{B- ^ D^K-) = A[B+ -^WK+) = Ab, 

A{B- Wk-) = Abz-, 

^ = Abz^, (3) 



where the complex numbers 



z± = rse 



i{SB ±7) 



(4) 



are the CP-violation parameters of interest, rs ^ 0.1 is 
the non-negative ratio between the magnitudes of the in- 
terfering b ucs and b cus amplitudes, and 5b is 
the CP-even phase difference between them. The magni- 
tude \Ab\ is measured [13] from the rate of the process 
B~ — > D'^K~, — > K~'!T~^, where contamination by 
the interfering decay chain B~ D'^K~, K~Tr~^ 
is doubly Cabibbo-supressed as well as rs-suppressed. 

We define the magnitudes Ap and RpAp to be the 
square roots of the total D° decay rates into F and F, 



Af = 0X00^^ = \/lXD°^^, 



F) 



1 



F). (5) 



The ratio Rp equals 1 for charge self-conjugate final 
states {F = F), but can in general have any non- negative 
value. Eqs. (5) ignore the possible impact of CP-violation 
in D decays. In addition, our use below of Ap and Rp 
will also ignore the effect of D° — D° mixing. It has been 
demonstrated [14] that these effects can be neglected for 
the purpose of measuring 7, as long as this is done con- 
sistently for the D mesons produced in the B decay as 
well as for those used to determine necessary D-decay 
quantities, discussed in Section III. Alternatively, pre- 
viously measured mixing and CP violation in D decays 
can be explicitly accounted for in the formalism [15]. For 
the purpose of the current discussion, it is sufficient to 
neglect these effects, as we do throughout this paper. 

We define the normalized amplitude distribution func- 
tions for the P-dependent charm meson decays, 



fEo{P) 



A{D° P(P)) 

a; ' 

A{do -> f (P)) 
ApRp ' 

A{D° F{P)) 
ApRp ' 

A{D° F{P)) 



(6) 



These functions satisfy the relations 



/^o(p) = /|o(p), f§AP)^fEo{P) (7) 

as a result of CP conservation in the charm meson decays, 
and are explicitly normalized, such that 



/ 



fEo{p)\ dp= / /^,(p) 



dP = l. 



(8) 



Accounting for the interference between the b ucs 
and b cus amplitudes in the B meson decays, the 
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amplitudes for the four full decay chains are obtained 
from Eqs. (3), (5), and (6), 

A{B-^F{P)K-) = Ao[fEo{P)+RFf^oiP)z-), 
A{B+ - F{P)K+) = Ao {Rpf^oiP) + fEo{P)z+) , 
A{B- ^ F{P)K-) = Ao (RpfUP) + f^,{P)z-) , 
A{B+ ^ F{P)K+) = Ao (/|o(P) + RFfEo{P)z+) , 

(9) 

where Ao = \Ab\Af- The observable P-dependent B- 
decay rates are the squares of these amplitudes, 



dr^iP) 



dP 



\A{B^ FiP)K^)\ 



d^±(P) _ l./z,± ^.mz^±^'2 



dP 



^ = \A{B^ ^ F{P)K^)\ 



(10) 



In the case F = F, namely, when the D decay fi- 
nal state is self conjugate, only two of the four equa- 
tions (10) are unique. These are the modes that have 
been studied experimentally so far [8-11, 16]. As men- 
tioned above, measurements of z± using F = K2.tt^tt~ 
and F = K^K~^K~ have been performed by analyz- 
ing only the P-dependence of the event distributions 
dr± (P) /dP, disregarding the total decay rates F^. Since 

fitting dT±{P)/dP in terms of rs, 7, and 6b leads to an 
average upward bias in rs when rs is of order its ex- 
perimental error, Refs. [8-10, 16] used the CP-violation 

parameters 



x± = 3?{z±}, y± = 9{z±}, 



(11) 



which are unbiased for this type of analysis. After these 
parameters are measured in the analysis of dr^{P) /dP, 
they are converted into (in general, non-Gaussian) confi- 
dence regions in terms of the "physical" parameters r^, 
7, and 6b- 

Here, however, we wish to focus on and generalize the 
approach used experimentally in Ref. [11] and first stud- 
ied theoretically in Ref. [17], by examining the additional 
information that can be extracted from the total decay 
rates and F^. The expressions for these rates are ob- 
tained by taking the squared absolute value of Eqs. (9) 
and integrating over all phase-space points, 

F^ = Al {1 + Rl\z-f - 2RFn{z*pZ-}) , 

Tl = Al {Rl + Iz+f - 2Rf^{zfz+}) , 

= Al {Rj, + \z^\^ ^ 2Rf^{zfz^}) , 

Rl\z+\^ - 2RF?li{z*pZ+}) , (12) 



At {I 



where 



ZF ^ - J fEo{P){fSo{P)ydP 

= - J /^o(P)(/|o(P))*rfP (13) 



is a measure of the interference between the and 
decay amplitudes into the final state F, averaged over 
the final-state phase space. The absolute value and ar- 
gument of 2;^ are, respectively, the coherence parameter 
and average strong phase of Ref. [17]. For the purpose 
of this discussion, it will be more useful to graphically 
think of 2:f as a coordinate-system offset parameter for 
z±. Methods to measure zf are outlined in Section III. 
The important point for now is that zf can be measured 
significantly more precisely than z± from high-statistics 
D decay samples, namely, 

CT2F<CT2±. (14) 

It is useful to represent z± in terms of the parameters 

1 



P± = z±- -5-ZF, 
p± = z±- RpZp. 



(15) 



We follow Ref. [11] in referring to p± and p± as the polar- 
coordinate parameters. This designation is motivated by 
the fact that measurement of the absolute decay rates is 
directly related to the radii \p±\ and |p±|, via the rela- 
tions 



^^(1 



Rl\p- 



- \ZF\-} , 

p+\- ^ RI\zf\'' 
Al (i?| + |p_p-i?||z^p 
Al{l + Rl\p^\^-\ZF\''). 



(16) 



Fig. 1 demonstrates the relationship between the polar 
coordinates p_, p_ and the Cartesian coordinates X-, j/_ 
for specific values of z_ and zp- The absolute values \p±\, 
\p±\ extracted from the total decay rates of Eqs. (16) 
yield two possible values for z_ and two for for a solu- 
tion of 7 with a four-fold ambiguity. In that sense, this is 
identical to measuring 7 with two, two-body D modes, as 
in the method of Ref. [18], whose discrete ambiguities are 
further discussed in Ref. [19]. In the case of multibody D 
modes, analysis of the distribution of events throughout 
the F phase space reduces the ambiguity to two-fold, in 
addition to improving the total precision [6]. In effect, 
the event phase-space-distribution analysis measures not 
only the absolute value but also the phase of p± and 
p± [11]. Combining the phase-space-distribution analy- 
sis with the total-rates analysis yields the most precise 
measurement of 7 for a given D decay mode. 



III. MEASURING zf 

A general approach for measuring the components of 
zp from decay rates of the ■0(3770) into neutral- Z? final 
states has been developed in Ref. [17]. Consider the case 
in which one of the ■0(3770) daughters decays into F{P) 
and the other decays into F'{P'), where the phase-space 
points P and P' do not have to be related. Due to the 
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FIG. 1: The relationship between the Cartesian and polar 
coordinates for zf ~ 0.5 — 0.2i, Z- = —0.1 + O.li, Rf ~ 
0.7. The solid (dotted) arrow corresponds to the complex 
number p_ (p_) of Eq. (15). A measurement of \p-\ (|p-|) 
implies that the true value of z- may be anywhere on the 
solid (dotted) circle in the {x-,y-) plane. The two crossing 
points of the circles are the possible solutions of z- . 



Next, we take F' to be a CP-even or CP-odd state, 
namely 



1 



=^ a/2 ^ ' 
The inverse relations of Eq. (21) yield 



(21) 



±1, 



(22) 



where the normaUzation condition Eq. (8) has been taken 
into account. Then Eq. (18) becomes 



fF,D± =2{1txf). 



(23) 



Eqs. (20) and (23) arc sufficient for obtaining xp and 
yp, the latter with a sign ambiguity. To resolve this 
ambiguity, we now take F' to be the 2-body state K~n~^. 
Eq. (18) then gives 

^F,K-n+ = ^F + ^lc-n+ 

- 2RFRK-n+ {xfXk--k+ +yFyK-w+){'^^) 



We have yet to determine xk- 
obtained from the rates 



T+ and yK-TT+- These are 



quantum numbers J^^ = 1 of the '0(3770), its two-D 

decay wave function is antisymmetric under exchange of 
the daughters, and hence must be {D^^D° - D°D°) /y^. 
The normalized event density in the PP' phase space is 



D _ A2 A2 
D^D'^ -"-F F' 



dT F,P'(PJ^') 

dPdP' 



fEoiP)RF'f§:{P') - RFf§oiP)fE'o{P'] 



,(17) 



K-TT+,K-n+ 



yK--!r+) ' 



^K-Tr+,D± = 2{1T XK--n+) ■ 



(25) 



as in Eqs. (20) and (23), respectively. Since K~n~^ is a 
two-body state, f^o '^^ and '^^ are numbers rather 
than functions. It then follows from Eq. (8), that ZK-Tt+ 
has unit magnitude, and the constraint 



(26) 



where Bj^„jj„ is the 0(3770) D^D" branching fraction. 
Integrating this expression over phase space yields the 
normalized rate. 



~ R _ 42 42 

^ no no^F F' 



F,F' 



'D«D 

Rp 



R 



R'p -\- R^ 



2RfRf>^{zfz*f,} 

2RfRf' {xfXf' + VfVf') j(18) 



where Eq. (13) was used, and we have separated zf and 
Zf' into their real and imaginary parts. 



ZF 
Zf' 



xf + iyF, 
xf' +iyF'- 



(19) 



By measuring the decay rate of Eq. (18) for different 
final states F', one obtains all the information about Zf- 
We begin with F' = F, for which Eq. (18) becomes 



tF,F = 2i?|(l-|zF| 



2Ri 1 



X F 



(20) 



resolves the ambiguity in yK--K+- Thus, it is possible 
to measure the real and imaginary parts of Zf with no 
ambiguities. 

Several studies [6, 20, 22] have shown that when ob- 
taining £)-decay parameters from '0(3770) decays, the 
expected error on the CP-violation parameters due to 
the finite ^'(3770) statistics is relatively small, given cur- 
rent CESR-c and _B-factory integrated luminosities. A 
detailed simulation study [22] has shown that in the 
phase-space-distribution analysis with F = Jf^Tr+Tr", 
the error on 7 due to the finite 0(3770) statistics is 
about four times smaller than the error due to the fi- 
nite B~ — > DK~ statistics in the currently available, 
1 ab"""^ i3-factory data sample. Measurements performed 
by CLEO-c with 818 pb^^ of e+e" ^(3770) data have 
yielded an estimated 7 error of 1° — 2° due to the mea- 
surement of the D KsTT^TT^ decay parameters [23]. 
CLEO-c has also measured zf for the modes X^tt+tt'^ 
'TT~^ , obtaining the preliminary results [24] 



\Zk- 
\zk- 



0.79 ± 0.08, arg{z;^- 
0.24 arg{zK- 



,+^0} = (197+^«) 



(16lJ)^ 
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Wc note that while these errors are large, their impact 
on the errors of \p±\, which arc the relevant CP-violation 
parameters for small- modes (see Section IV B) is sup- 
pressed by Rp, as seen from Eqs. (15). 

Furthermore, the newly launched BEPC-II charm fac- 
tory, with a design luminosity almost twenty times that of 
CESR-c, will be able to supply the charm data needed to 
match the large B samples that will be collected at LHCb 
and possibly at a proposed e+e~ "super B factory". We 
conclude that the error on z± in a decay-rate analysis of 
the type presented in this paper will be dominated by 
the experimental error on \p±\ and not by knowledge of 

ZF- 



A. Self-Conjugate Modes 

So far, all the multibody iP-decay modes studied exper- 
imentally within the context of B^ — > DK^ have been 
charge self conjugate, i.e., F = F. From Eq. (5), one sees 
that self-conjugate modes satisfy Rp = 1. In addition, 
Eq. (7), together with the condition F = F, implies 



.fUP) = .f^o{P)- 

As a result, such states satisfy 
2/F = 0, 



(27) 



(28) 



as we demonstrate by dividing the phase space of F 
into two equal- volume regions V and V, such that ev- 
ery point P G V is related to a point P G V hy the 
CP transformation. For example, in a three-body de- 
cay of the type D a'^a~b°, the division is along the 
line {pa+ + Pb")"^ = {Pa- +P6°)^) where pj is the four- 
momentum of particle j. Such a division can be per- 
formed for any multibody final state that is self conju- 
gate, regardless of its particle multiplicity. Then 

VF = ^y^fUP) (/io(P))*dP 

jEoiP) {f^oinYdP 

(29) 



Using Eq. (27), the second integral in Eq. (29) can be 
written as 



/|o(p) {fEo{p)ydP 



(30) 



The integrand in Eq. (30) is the complex conjugate of the 
integrand of the first term in Eq. (29). Therefore, their 
imaginary parts cancel in the sum, yielding yp = 0. 

In the 7-relatcd measurements performed so far with 
B"^ DK^ and D F decays into a multibody, self- 
conjugate state, a particular model was assumed for the 
functional form of f^o{P). The parameters of the model 
were obtained by fitting the phase-space distribution of 



£)" F decays, where the flavor of the D° was tagged 
by its production in the decay _D*+ Z?°7r+. In this 
case, Eq. (27) guarantees that zp can be fully determined 
by inserting the model f^o{P) into Eq. (13), as done in 
Ref. [11]. The same cannot be done for modes that are 
not self conjugate, where one must resort to the use of 
V'(3770) decays. 



IV. EXPERIMENTAL SENSITIVITIES 



Due to the linear relationship (16) between the experi- 
mentally observable rates and |p±p, |p±P, these squared 
radii are the unbiased CP- violation parameters of choice 
for the rates analysis, given that decay rates can almost 
always be obtained from reasonably unbiased estimators. 
If the errors on \p±\^ and |/0±P are significantly smaller 
than the values of these parameters, then their roots |p±| 
and \p±\ are also unbiased parameters. 

In terms of the errors on the rates, the errors on \p±\ 
and \p±\ are 



(TpF 



Hp- I 



Hp- I 



^|p+l 



^ + 



21^' 
2\p+\AlRy 



(31) 



where we have used Eq. (14) to neglect the error on \zp\. 
As a result of Eqs. (14) and (15), the errors on \p±\ and 
\p± \ are similar in magnitude to the errors on x± and y±. 
Therefore, studying Eqs. (31) provides a simple means to 
compare the 7 sensitivity of a rates analysis using any fi- 
nal state F to the sensitivity of the current-best measure- 
ment, namely, that of x± and y± from the phase-spacc- 
distribution analysis oi F = K'^tt'^tt~. In what follows, 
we make quantitative estimates of the errors on \p±\ and 
|P±|- 



A. Self-conjugate modes 

As a result of Eq. (28), Eq. (15) simplifies to 

P± = P± = z± - xp (32) 

for self-conjugate modes. Therefore, the two circles of 
Fig. 1 collapse onto each other, and the rates measure- 
ment of the radii \p±\ is no longer sufficient for fully de- 
termining z±. This is hardly a problem, for two rea- 
sons. First, the phases of p± may be determined from 
the event-distribution analysis, as was done in Ref. [11], 
yielding a measurement of z± whose precision is enhanced 
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due to the use of all available experimental information. 
Second, as stated in the introduction, precise knowledge 
of 7 can in any case be obtained only by combining many 
measurements of parameters related to 7. Therefore, 
measurement of |yO±| helps reduce the overall error on 
7, even if it is not sufficient for extracting 7 without in- 
formation obtained from other 7-related measurements. 

Since xp is well known, it is useful to estimate the 
errors on \p±\ for relevant D decay modes, as they will 
correspond closely to the errors on z±. Since we are deal- 
ing with the case Rp = I, Eqs. (31) become 



TABLE I: Values of the inputs to Eq. (33) and the expected 

errors (7\p^ \ for different D decay modes and a B-factory data 

sample of lO'^ e^c^ BB events, calculated with the CP- 
violation param(>t(>rs (jf Etj. (-i!)). 



Mode 


Xf 


r_ 


r+ 


\P-\ 


\p+\ 






7r"'"7r~7r" 


0.85 


0.13 


0.10 


0.75 


0.94 


0.07 


0.06 




-0.02 


0.05 


0.05 


0.13 


0.08 


0.19 


0.32 


KlK+K- 


-0.31 


0.09 


0.10 


0.41 


0.23 


0.12 


0.20 


K+K-n° 


0.20 


0.11 


0.11 


0.12 


0.29 


0.48 


0.19 



1 



\P±\ 



2|P±| 



(33) 



where we have used T± = FfJ. = FfJ., and the second 
equality of Eq. (33), obtained from Eq. (16), conveniently 
relates the \p±\ errors to the relative errors on the signal 
branching fractions. We rely on previous "reference" ex- 
perimental studies of the relevant decay modes to obtain 
these relative errors for any hypothetical value of \p±\- 
Suppose that in a reference measurement performed with 
i3-factory data of integrated luminosity L, one observed 
N± signal FK^ events, from which the rates 

were determined and the CP-violation parameter values 
were calculated. Let N = N++N^. Then the num- 
bers of signal events that would be observed in an experi- 
mentally identical, hypothetical measurement of luminos- 
ity L given hypothetical values |p±P for the CP- violation 
parameters, are 



N± = Nr± — , 
L 



where 



l + \p±\'-x% 



2x% 



(34) 



(35) 



^ given the 



is the ratio between the value of the rate F;'^ 
hypothetic;al parameter values p± and the sum of the 
rates + t'^ measured in the reference measurement. 
The second equality in Eq. (35) arises from Eqs. (16). 

We assume that the error on the number of events N 
in the reference measurement can be written as the sum 
in quadrature of a Poisson signal part and a background 
part, namely. 



4 = ^ + 4,6,d- 



(36) 



Using this relation and the published reference- 
measurement quantities N and aff, we obtain the back- 
ground contribution to the error, which we assume to 
be CP symmetric. Then the errors on the numbers of 
events in the hypothetical measurement, in which N± 
will be observed, are 



(37) 



where the statistical assumption leading to Eq. (36) was 
again used. From Eqs. (34) and (37), we obtain the rel- 
ative branching-fraction errors for the hypothetical mea- 
surement, 



0T± 



TV 



Nf± 



Nf4 



(38) 



In Table I we report xp, \p±\, and cr|p^| for sev- 
eral three-body D decay modes, assuming a data sam- 
ple of 10® e~^e^ BB events, similar to the cur- 
rently available i3-factory sample. We obtain the val- 
ues of xp from Eq. (13), using the Dalitz-plot distribu- 
tions f^o (P) , whose parameterizations are reported in 
Refs. [11], [8], [10], and [25] for the D-decay final states 
tt+tt-ttO, K^n+TT-, K^K+K-, and K+R-tt^, respec- 
tively. We also obtain N, ct^, and L from these refer- 
ences, except for K~^K~Tr°, where we estimate N and aj^ 
from their values in 7r+7r~7r'^ [11], taking into account the 
ratio of branching fractions B{D" ^ K+ K-tt^)/B{D^ -f 
7r+7r~7r°) [26] and an assessment that the background 
yield in K+K'-k^ will be 20% of that in -k+'k'-k^. Since 
extraction of |/5± |^ from the total rates has been reported 
only for the F = tt+tt^tt'^ mode [11], we take |p±| = |p±| 
when evaluating f± for all other modes, for lack of a 
better value. We take the hypothetical CP-violation pa- 
rameter values \p± I from the averages of the values of x± 
and y± reported in Refs. [10] and [16], 

x- = 0.097 ±0.034, 

y- = 0.054 ±0.058, 

x+ = -0.087 ±0.031, 

y+ = -0.038 ±0.042. (39) 



The errors of Eq. (39) reflect the sensitivity of a mea- 
surement conducted with 1.04 x 10® e+e~ BB events, 
comparable to the value used to produce Table L 

One can see from Eq. (33), that the error a'\p_^\ is small 
when jp-tj is large. Large requires xp to be large, by 
virtue of Eq. (32) and the smallncss of demonstrated 
in Eq.(39). We note that some insight into the value of 
Xp for a particular mode can be obtained by studying 
the distribution of events in the D*'-decay Dalitz plot, 
since generally, high level of apparent symmetry under 
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the exchange of the two charged particles leads to a high 
value of xp. 

It is evident from Table I that of the three-body 
modes studied here, only Xc^^+tt-ttO is large enough for 
Eq. (33) to yield \p±\ errors that are competitive with 
the errors of Eq. (39). In particular, the high-statistics, 
low-background mode -ftTgTr+Tr" ends up having large \p± \ 
errors due to the very small value of xj^o^+^-. On the 
other hand, we expect that the methods for suppression 
of the significant background in 7r+7r~7r°, which were first 
developed in Ref. [31], will improve in upcoming analy- 
ses 

simple extrapolation shown in Table I. 



That should reduce for this mode below the 



B. Non-self-conjugate modes 

We proceed to estimate the errors on p± and p± in TP- 
decay final states that are not self-conjugate, i.e., F ^ F. 
As in the procedure leading up to Eq. (33), we replace 
Al in Eqs. (31) using Eq. (16): 



(TpF 

2\p_\Rl 



^\p-\ = 



2|P+| 



a. 



^\p+\ = 



2\p. 



(40) 



where 

l^p^ = l + i?||p±p- |ZFP, 

z>p^ = R% + \p±\^ ~ RI\zf\\ (41) 

As in Eq. (38), the relative errors in Eqs. (40) are ob- 
tained from the number of signal events , and 
their errors, observed in existing reference measurements. 



CTpF 



(42) 



with an analogous expression for F, where by analogy 
with Eq. (35), 







^P- + % 








= 




y-p^ + v-p_ 




Vp+ 


Vp^ + Vp_ 



(43) 



As in Section IV A, the symbols p±, in Eq. (43) re- 
fer to the CP-violation parameters extracted from the 
reference measurements and . If the total rates 
were not used to extract CP-violation parameters, one 
can naively take p± and from Eq. (39) for the pur- 
pose of performing this error estimate. 

Let us consider this error estimate in the case of the 
non-self-conjugate, three-body final state F = KgK~'K~^. 
With as little as 5% of their currently available data sam- 
ple, the BaBar collaboration has performed a preliminary 
analysis of this mode's Dalitz-plot amplitude-distribution 
functions f^o{P) and .f^a{P) [27], from which we com- 
pute \zf\ = 0.47. The ratio Rk°k-t!+ = 0-68 is easily 
extracted from the results reported in Ref. [27]. With 
Rk'^k-it+ being different from 1 yet of order 1, this mode 
is in a class of Cabibbo-suppressed decays expected to ex- 
hibit large interference between the b — > ucs and b — > cus 
decays [28]. Unfortunately, as we show below, the combi- 
nation of a small branching fraction and a medium-sized 
\zf\ render KgK~7T^ unattractive for extracting 7 via 
the total-rate method. 

In addition to i?x°K-7r+ and \zkok--k+\j calculation 
of all four errors of Eq. (40) also requires knowledge 
of arg{2^o^-^+}, which has not been measured. How- 
ever, a rough estimate of the CP-parameter errors shows 
them to be comparable to those of the K'^K~^K~ mode, 
due to the following two observations. First, the com- 
bined branching fraction B{D° K°K-Tr+) + B{D° 
K°K+Tr-) is approximately 85% of B{K°K+K-). One 
therefore expects the relative error on 
somewhat larger than that on N^°s^*^ . Experimental 
details, such as kaon vs. pion multiplicities and combi- 
natoric background under the larger K^Kn Dalitz plot, 
slightly increase our expectation for the ratio between 
the relative errors on N^s^^ and N^s^ ^ . The sec- 
ond observation is that l-s^^o^-jr+l is about 50% larger 
than Xx°x+K-- Combining these two competing ef- 
fects, we conclude that the errors on components of 
the CP-violation parameters obtained from K^Kit and 
K^K^K~ should be of similar magnitudes. As seen in 
Table I, this implies error values that are too large to be 
of practical interest. 

We note that Eqs. (40) also hold for Cabibbo-allowed 
final states involving a single charged kaon, such as F = 
K~-K^-K^, for which Rj^-^+^a w 0.05 [29] (where we have 
ignored the effect of D'^ — mixing [14]). Eqs. (16) 

show that in this case, the sensitivity of and F^ to 
the CP-violation parameters is suppressed by i?|-_^+^o, 
making these rates \iseful for obtaining Ab , as mentioned 
in Section II for the K~tt^ decay. However, the 

absolute rates and F^ do provide a good measurement 
of \p±\. Searching for these decays in a data sample of 
226 X 10*5 e+e" BB events, BaBAR [30] has put an 
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upper limit on the ratio 



Rads 



■'P+ 



(44) 



p+ 



for which the central value obtained was Rads — 
0.013 The rates that appear in the numerator 

of Eq. (44), to which we refer as the ADS rates [18], are 
suppressed by factors of second order in the small param- 
eters rs, Rx-T!+T!° relative to the rates in the denomina- 
tor. The error on Rads is dominated by the statistical 
errors on the ADS rates. To properly account for this 
when calculating the relative errors on the ADS rates, 
we evaluate Eq. (42) with 



(45) 



instead of the expressions in Eq.(43), and take N^'^'^ 
to be the number of ADS events detected in Ref. [30], 
namely, 19±10, where the 10-event error is obtained from 
the naive average of the positive and negative errors on 
Rads- 

The resulting errors on |p±| are shown in Fig. 2, cal- 
culated with Eqs. (40) for different values oi zp- As in 
the case of Table I, wc have assumed a data sample of 
10^ e~^e~ BB events and the CP-violation parame- 
ter values of Eq. (39). The errors reach values as low as 
0.016 and as high as 0.035 (0.045) for (|p+|). We 
sec that at least one of the errors is smaller than about 
0.025 for any value of zp- For the CLEO-c central values 
of Zk-t,+-kO [24], we find a\p^\ « 0.02. 

These results suggest that one can expect measurement 
of the CP-violation parameters with F = K~n^TT^ to 
yield errors that are very competitive with the current- 
best measurement, Eq. (39), once the luminosity is high 
enough for observation of the ADS decays. 




FIG. 2: The errors on |p_| (top) and jp+j (bottom) as func- 
tions of the absolute value and phase of zf for F = K~n~^n'^, 
calculated for 10^ e+e" — > BB events with the CP-violation 
parameters of Eq. (39). 



V. DISCUSSION 



Of the self-conjuage final states studied quantitatively 
here, the errors obtained from 7r+7r~7r° are the smallest, 
due to the large value of xp = ^{zp} = 0.85 in this 
mode. The errors are expected to decrease beyond the 
estimate shown in Table I, as background suppression 
improves in subsequent analyses of this mode. By con- 
trast, the final state K'^tt^tt^ , which thanks in part to its 
large branching fraction and high purity has yielded the 
most precise phase-space-distribution measurements of 7 
to date, has a very small xp, rendering its absolute decay 
rates poor measures of the CP-violation parameters. 

Our calculations show that measuring \p±\ with the 
final state iir~7r+7r° can yield very small errors, smaller 
than or of similar magnitude to the errors from the phase- 
space-distribution analysis of K'^tt~^tt~. We note that 
similar precision may be obtained with the two- and four- 
body final states K~tt~^ and K~TT'^Tr~Tr~^ , whose study 
is outside the scope of this paper. 



The results presented here cover the major three- 
body D-decay final states with known and significant 

branching fractions. It is possible that the absolute de- 
cay rates into some of the higher-multiplicity states will 
also turn out to yield competitive errors on 7. Among 
the Cabibbo-suppresscd modes, this includes the final 
state K^TT^TT^ , whose phasc-spacc-distribution anal- 
ysis has been studied in simulation [12], and 27r^27r~. 
The Cabibbo- favored mode — > K'^tt'^ 'k~ has a large 
branching fraction, (5.3 ± 0.6)% [29], and may there- 
fore be attractive for both phase-space-distribution and 
absolute-decay-rate analyses. Since almost half the rate 
is due to the resonant contribution ii'*~(892)p+, the 
phase-space distribution is highly asymmetric under ex- 
change of the two charged pions. Therefore, it is unlikely 
that Xp is large for this mode. Nonetheless, given the 
large branching fraction, even xp as small as 0.1 could 
make this mode attractive for studying 7. 
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VI. CONCLUSIONS 

We have studied the use of the absolute DK^ 
decay rates, where the D decays to a multibody final 
state, for obtaining information with which to improve 
the overall knowledge of the CKM unitarity-triangle 
phase 7. This information is complementary to that 
obtained from other 7-related mcasiircments, including 
analysis of the _D" — D° interference pattern seen in the 
phase-space distributions of the D decay products. We 
have developed a formalism for estimating the error on 
the CP- violating parameters \p± \ and |jO±|. The parame- 
ter that most strongly affects the errors is zp of Eq. (13). 
We have evaluated zp for three-body D final states for 



which the necessary input information is available, and 
have estimated the errors on the CP-violation param- 
eters for these self-conjugate modes and for the modes 
K'^^K--K+ and K-tt+-k°. 
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